Indications of coherence-incoherence crossover in layered transport 
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For many layered metals the temperature dependence of the interlayer resistance has a different 
behavior than the intralayer resistance. In order to better understand interlayer transport we 
consider a concrete model which exhibits this behavior. A small polaron model is used to illustrate 
how the interlayer transport is related to the coherence of quasi-particles within the layers. Explicit 
results are given for the electron spectral function, interlayer optical conductivity and the interlayer 
magnetoresistance. All these quantities have two contributions: one coherent (dominant at low 
temperatures) and one incoherent (dominant at high temperatures). 

PACS numbers: 71.38.-k,71.38.Ht,72.90.+y 



Many of the most interesting strongly correlated elec- 
tron materials have a layered crystal structure and highly 
anisotropic electronic properties. Examples include the 
cuprates [jj], g, |), colossal magneto-resistance materi- 
als 0, U , organic molecular crystals [|J [?J , and strontium 
ruthenate H . One particularly poorly understood prop- 
erty of these materials is that often the resistivity perpen- 
dicular to the layers has quite a distinct temperature de- 
pendence to that parallel to the layers |J . This is in con- 
trast to what is expected for an anisotropic Fermi liquid: 
the parallel and perpendicular resistivity then have the 
same temperature dependence, being determined by the 
intralayer scattering rate, r(T). (This result holds even 
when r(T) is larger than the interlayer bandwidth Jlo[). 
In many of these materials the interlayer resistivity is a 
non-monotonic function of temperature with a maximum 
at some temperature Tj nax . In some of the materials the 
intralayer resistivity also has a maximum as a function of 
temperature, but at a higher temperature T|| nax > TJ nax . 
Previously it has been suggested that the maximum in 
the temperature dependence of the interlayer resistivity 
is associated with a crossover from coherent interlayer 
transport at low temperatures to incoherent transport at 
high temperatures. There is no consensus as to what 
actually determines T™ ax . Recent angle resolved pho- 
toemission spectroscopy (ARPES) experiments on two 
different layered cobalt oxide compounds JllJ found that 
one only observed peaks in the electronic spectral func- 
tion (corresponding to coherent quasi-particle excitations 
within the layers) below a temperature T coh that was 
comparable to T™ ax . 

Although many theoretical papers have considered the 
problem of incoherent interlayer transport (see for exam- 
ple Ref. [jj] and the references given in Ref. [l(]), we are 
unaware on any theory which starts with a many-body 
Hamiltonian and produces the three temperature scales, 
T coh , J"™ ax ; and T| nax . In this paper we consider a simple 

microscopic model which exhibits T coh — T™ x < T™ ax . 
The model is a layered version of Holstein's molecu- 
lar crystal model where the electrons strongly couple to 
bosonic excitations to produce small polarons. It should 



be stressed that we are not claiming that the charge 
transport involves small polarons in all of the above ma- 
terials . Rather, we suggest that this model can provide 
insight into the relevant physics associated with these 
temperature scales. 



We start with a Holstein model 13 for an infinite lay- 
ered system where the electrons interact with dispersion- 
less bosons. The Hamiltonian is 



h = ^ e ° c l c i + ^0 X! a i a i + X! t "i c l c i 

i i 

+Mj24c t (a l + al), 



<ir/> 



where e° is the on-site energy, ujq is the characteristic 
frequency of the bosons, U v is the hopping integral be- 
tween nearest-neighbor sites i and 77, M is the coupling 
between the bosons and the electrons. We introduce a 

dimensionless coupling g = {^j^j • We have that g > 1 

in order for small polaronic effects to be important. Since 
we want to study layered systems we split the hopping 
into parallel and perpendicular to the layers, tu and t± 
respectively where t\\ t±. This enables us to write the 
Hamiltonian in a way more adapted for the layered case, 
shown in Fig. ^. The nature of the transport depends 
on how iii and t± compares with T, the scattering rate 
due to the bosons. We assume that that T > t±, so that 
the inter-layer transport can be described by considering 
two decoupled layers. The Hamiltonian can be specified 
for this system. Two layers are coupled with a hopping 
Hamiltonian. Within each layer the electrons can hop 
but there is a coupling to a bosonic degree of freedom in 
each layer with characteristic frequency Wq. The bosons 
can be phonons, magnons, plasmons, or particle-hole ex- 
citations. We only consider a single frequency loq for 
reasons of simplicity; it allows us to express some of our 
results in a analytical form. 
We then use the Hamiltonian: 

H = H-i + Ti-2 + T~Lt 
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where 
Hi 



e°c\ci + tkJo at Of + 1\\ ^ c* 
+M^cfci(ai + at), 



<J7)> 



W 2 = 



X e °4 d j + huj o Y, b h + h Y, d Uj 

j j <jS> 



n t = t ± j2( c l d *+ h - c -)- 

i 

Here, and below, c, i, a, 1 refers to one layer, and d, j, b, 2 
to the other one. First we focus on the properties of the 
two separate layers. We perform a Lang-Firsov transfor- 
mation jlij to remove the coupling of the electrons to the 
bosons. Then c, — > & 



Xi = exp 



and Oj — > ai — —c\ci. The Hamiltonian is transformed to 

= e s He- s where S = ^ ^ 4 c *( a ! ~ Thi s diag- 
onalizes the electron-boson part of the Hamiltonian, but 
introduces extra X-operators ]l5| ] in the hopping parts of 
the Hamiltonian. For example: 



M < h 



(i) 



W t =t x ^(ctxt^. + hx.), 



(2) 



<ij> 



The intralayer hopping terms can be treated by adding 
and subtracting to the Hamiltonian a term 



(3) 



<y> 



where (..) denotes a thermal average over boson states 
and this term describes a tight-binding band of small 
polarons within each layer 114, [t6) 

e k = e° - e-^ 1+2nB h\\ [cos(k x a) + cos(fc y a)], (4) 

where a is the lattice constant within the layers and 
ns{T) = (exp(hu!o/kBT) — is the Bose function. 




FIG. 1: The two layers with polarons are coupled by a hop- 
ping Hamiltonian. 



There is then a residual interaction |16| between the 
polarons and the bosons which is described by 



n 



p—b 



(5) 



and leads to scattering of the small polarons. The first 
non-zero contribution to the imaginary part of the po- 
laron self-energy, E, comes when the polaron emits one 
boson and absorbs one boson. If the energy dependence 
of the the density of states (DOS) is neglected one finds 
that liUl. 



T(T) = Wg 2 n B (T)(l + n B (T)), 



(6) 



where W is the renormalizcd bandwidth, W = 
4<||e-9( 1+2 " B ) Qjnjffil. 

Note that the small polarons are composite particles 
(quasi-particles) . They consist of an electron bound to 
a "cloud" of bosons. This coherent quantum state can 
move freely within the layers producing coherent charge 
transport. In contrast, ARPES involves ejection of elec- 
trons rather than polarons from the crystal. Similarly, 
the interlayer charge transport involves the tunneling of 
electrons between layers. In order for this to occur the 
bosons bound to the electron in the polaron must be re- 
moved, the electron tunnels, and a new set of bosons are 
bound to the electron. 

Greens function within a single layer. Let us start by 
calculating the Matsubara Green function (GF) within 
one layer, ignoring the coupling between the layers. After 
the Lang-Firsov transformation we calculate the small 
polaron GF 

G°(k,r) = -iQ( T )±Y^ <Ri '~ Ri) ( T T~ c i'(T)4m 



= -iQ{T)e ( ^- ir)lr/h , 

where 0(r) is the step function, T T is the time ordering 
operator, is the position of lattice site i. T is the 
inverse lifetime associated with the interaction, Eq.(^|), 
and there are N lattice sites in each layer. 

The electron GF, G(k, iuj n ) involves a convolution of 
the polaron GF with the Fourier transformed (denoted 
J-) A-operators [[L7| 

G(k^„) = l Y ^{<^MA J t (0))}G°(k',c„0e 4(k - k ' HR ''- R 

u) n i ,k' 

Performing the summation over RLj/ — R.^ we obtain 

1 



G(k,iw„) = e -f( 1 + 2 "«) 



iuj n — ek + *r 

- 1 



iui n — ef + iT 



k',2/0 



2g\Jn B (l + n B ) 



-lHuoP/2 



iu> n ~ £k' + IhujQ + iT 



(7) 
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FIG. 2: Energy dependence of the electron spectral function 
at a wave- vector on the Fermi surface, using a constant DOS. 
The product of the spectral function A(k.F, e) with the Fermi- 
Dirac distribution function /(e) is shown because this can be 
compared with ARPES spectra. Note that the well-defined 
quasi-particle peak which occurs for ksT hu>o disappears 
at higher temperature. The results are shown for an electron- 
boson coupling of g = 1 . The inset shows the same thing for 
a smaller bandwidth. 



Note that the GF is a sum of a coherent and a inco- 
herent part, i.e., the terms on the second and third lines 
are independent of k. Eq. (Q) is a generalization to non- 
zero temperatures of the GF found by Alexandrov and 
Ranninger p7[ . 

In Fig. pPwe plot the electron spectral function, 
A(k,uj) = Im[G(k,w)], resulting from the GF, Eq.(Q), 
for different temperatures. With increasing temperature 
the boson modes become populated, ns(T) increases and 
the spectral weight shifts from the coherent part of the 
spectral function to the incoherent part. 

Qualitatively similar behavior was seen in recent 



ARPES [pi measurements. This behavior does not 
change much qualitatively when g is changed. From plots 
we estimated that the crossover takes place at 



coli 



(8) 



This can also be justified using Eq.(Q) when W < frojo. 

Interlayer conductivity. Let us derive an expression for 
the current perpendicular to the planes. We assume that 
the hopping between the layers only takes place between 
nearest neighbors (see Fig. |lj). At applied field V the 
current is given by 



I^{eV) = ^Im 



dr 



ieVr 



(TVi(r)jt(O)) 



(9) 



where j is the current operator, fj, and v are directions 
in the crystal. 

Using the Hamiltonian above and Eq.(|J) we have, to 
second order in t±, 



I±(eV) 



2e. 
~h 



dre 



ieVr 



^ 2 E/ 

x {T T c](r)d n (r)dl(0)c n (0)) (t t Y£ {r)Y 3 (r)Xj(0)X jl (0)) , 



where d is the distance between the 2 layers. 

We assume that the bosons in separate layers are in- 
dependent of one another. Hence, in Eq.([lO|), we can de- 
couple the X and Y polaron operators corresponding to 
the first and second layers. This means that the Fourier 
transformed averages of the electron operators gives rise 
to two GFs. These GFs describe polarons bands within 
each layer. 

After some algebra we can write the current as 



I ± (eV) = |^ e - 29(1+2nB) Jj~j££ ^ ?(kj e)A o {K e + eV) [f{e) /(£ + eV)] 



+ 



- 1 



(l 4 9v /n B (l + n B ) 

CO 

X! 7 ' ["WM 1 +n B ) 



dc_ 
2^ 



-lhui /3/2 



J2 ^?(k, e) ^(P, e + eV) [/(e) ~ f(e + eV)] 



<k_ 

2tt 



^?(k, e) £ A° 2 (p, e + eV + lhu ) [/(e) - /(e + eV + lhu; a 



(10) 



p belongs to the first (1) layer and k to the second (2) laron GFs in each layer respectively, and J; is a modified 
layer. A® and A\ are the spectral functions for the po- Bessel function of order I. 
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Note the similarity in structure between Eq.([7|) and 
Eq. (|lOj) . The expression for the current, Eg. (plf) , has 
a contribution from coherent and one from incoherent 
transport. The first line corresponds to tunneling where 
the momentum of the polaron parallel to the layers is 
conserved. In the second and third terms the intralayer 
momentum is not conserved. The third line has a dif- 
ference energy of IHloq between the polarons in the two 
layers because there is a non-zero difference between the 
net number of bosons that are absorbed and emitted in 
the tunneling event. At low temperature the coherent 
part dominates but at high temperature (fc^T > Hloq) 
the incoherent mechanism of transport will dominate. 



Thus, there is a crossover from coherent to incoherent 
transport. 

Then, the inter-layer conductivity is easily obtained 
This should be multiplied with the 



d(eV) 



eV=0 



number of sites in one layer. The interlayer resistivity, 
measured in experiments, is simply crj 1 . 

By using the fact that when T <C W, J2k e) = 
-D(e), where D(e) is the DOS, and assuming that terms 
containing the derivative of the spectral function is small 
dA{e + eV)/dV < A(e + eV)D{e) We can write the con- 
ductivity as: 



J 



ax = ^-t\d 2 e- 2 ^+^B) 



* W ' 2 de , , 

-W/2 27T 



igy/n B (l + 



h [4gy/n B (l + n B )\ e -«/2 / ^D(e)D(e + lfiw ) 



riB) 



-dm 

de 



de „, , „, ,, , —df(e + Itkoo) 

de 



(11) 



In general, the interlayer conductivity for decoupled 
layers is [ju| , 



It is interesting to note that a 



coll 



o- mcoh actually 



4f 
de 



(12) 



where A(\t, e) is the electron spectral function for a single 
layer. Directly substituting Eq.(^) in this we can obtain 
Eq.©. 

We can estimate the temperature of the crossover in 
the conductivity by having equal contribution from the 
coherent and the incoherent parts. If we look at the con- 
ductivity there is a minima (maxima in the resistivity), 
corresponding to the crossover (see inset in Fig. ||). Ig- 
noring the contribution from the I ^ terms in Eq.([lO|) 
we can get an approximate expression for the crossover 
temperature for the interlayer resistivity: 



„ 0.59^. (13) 

<f¥ g g 



k B T'l iax 



This expression compares quite well to the crossover tem- 
perature extracted from a numerical plot of the resistivity 
versus temperature in Fig. ^. Hence, we see that T coh and 
jimax are com p ara bi ep At these temperatures, T ~ 0.3PF 
A more detailed analysis shows that the I ^ terms in 
Eq.(|l0|) contributes linearly in W. So the estimate for 
the crossover, Eq.(|l3|), should only hold when W < hw , 
but extensive numerical work showed that it has a wide 
range of validity. 



occurs at a temperature lower than T™ ax , as can be seen 
in the inset in Fig. |[ This is because an almost constant 
contribution from the I ^ terms. 

Intralayer conductivity. When we calculated the cur- 
rent between the layers we assumed that the small po- 
laron band was well-defined. However, it is known that 
even within the layer, as the temperature increases, there 
is a crossover from coherent to incoherent transport, oc- 
curring at Ty 118 *, and there is a maximum in the resis- 
tivity associated with this crossover [fl3|| . An important 
question is the size of T™ a * relative to T™ ax . 

We are unaware of a systematic method of calculating 
the full temperature dependence of the intralayer con- 
ductivity. At low temperatures it can be calculated by 
considering the band transport of local polarons. At high 
temperatures one can consider the diffusion (or hopping) 
of localized polarons. Extrapolating the results to inter- 
mediate temperatures provides a means to estimate the 
crossover temperature and the regions of validity of the 
different expressions. 

At low temperatures (r <C W) the transport in the 
layers are coherent. We assume that we have well devel- 
oped quasi-particles in the layers and the bosons acts as 
a small perturbation. The conductivity is given by: 



e 

2^2 



de 



,2, «( k ) 2 r/ n f d f 



We use a tight binding model for the electrons to calcu- 



5 



late the velocities, and the decay T is taken from Eq. (g) 
above. 

At high temperatures (T ^ W) the electron is trapped 
by the formation of the polaron. The electrons are local- 
ized at the lattice sites and the concept of a wave-vector 
for the small polaron is meaningless. The intralayer hop- 
ping term in the Hamiltonian should then be treated as 
the perturbation. The intralayer conductivity at high 
temperatures from Eq.(^) depends on the polaron GF, 
solved when the hopping term is the perturbation. How- 
ever, due to the lack of dispersion of the bosons the time 
integral in the conductivity diverges Jl5| . We overcome 
this problem by including a finite self energy in the elec- 
tron GF. To evaluate this we sum a series of the lowest 
order diagrams. We have to take special care about the 
number of paths on the 2D lattice. If G(u) is the local 
electron GF it has a self energy that satisfies 



At\G 



1 + t\G 2 



(1 



f|G 2 ) 3 



(14) 



This is a self-consistent equation for the real and imagi- 
nary part of the self energy that have to be solved simul- 
taneously to self consistency. This was done on a square 
lattice where the polaron can hop to four different sites 
each time. The conductivity for the low and high tem- 
perature regions were plotted and the maxima in pn was 
extracted. 

There is a crossover also for the intralayer resistivity. 
From Fig. [|, the crossover occurs roughly when 



hljjn 



(15) 



using this in Eq.(||), means that the crossover occurs 
when r ~ 8W. The intralayer crossover occurs at higher 
temperatures than the interlayer crossover, therefore the 
assumption made above that for the interlayer calcula- 
tion we have well developed quasi-particles within each 
layer is justified. 

Optical conductivity. The frequency dependence of 
the interlayer optical conductivity, erj_(w), has been sug- 
gested to be a probe of interlayer coherence in the metal- 
lic state Jil| . The optical conductivity can be found from 
Eq.(^|), by letting eV — > uj, and calculating the derivative. 

Fig. |i] shows how at low temperatures there is a well- 
defined Drude peak at zero frequency due to coherent 
interlayer transport of small polarons. The width of 
this feature is approximately T. Note that this feature 
occurs even though T > t±, as has been pointed out 
previously Q. As the temperature increases the spec- 
tral weight of this feature decreases and is replaced with 
a broader feature associated with incoherent interlayer 
transport and with a width that is determined by the 
small polaron bandwidth within the layers. The incoher- 
ent part becomes narrower with increasing temperature 
because of the polaron narrowing of the bands. Changing 
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FIG. 3: Peak temperature in the resistivity as a function of 
coupling. Results obtained using a tight binding DOS (no 
major change was seen when using a constant DOS). The 
two sets of data-points are for the intralayer and the in- 
terlayer crossover temperatures respectively. The intralayer 
crossover occurs at much higher temperatures than the in- 
terlayer. The inset shows the interlayer conductivity as a 
function of temperature when g — 1 and W° = 77/kJo- 
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FIG. 4: Optical conductivity divided into the two contribu- 
tions, coherent and incoherent, plotted for two different tem- 
peratures. In the lower left panel we plot the coherent part, 

and in the right the incoherent part, <tq — 2e 2 d 2 I £J- ) jh. 



g and i|| does not qualitatively change this temperature- 
dependent crossover. 

Magnetoresistance. If we apply a magnetic field, B, 
parallel to the layers (the x-y plane) we have a orbital 
effect on the paths of the electrons. This can be de- 
scribed by a shift in the Bloch wave vector, k — > k— f A, 
where A is the vector potential for the magnetic field. 
For a magnetic field in the x direction, when an elec- 



6 



tron tunnels between adjacent layers it undergoes a shift 
in the y-component of its wave vector by -dB JhJ. In 
the general expression Eq. (|l2|) A 2 (k, e) is replaced with 
A(k, e)A(k + j^dBy). However, since the incoherent part 
of the conductivity contains a summation over k-space 
and is independent of k, this will be unaffected by the 
magnetic field. Thus, we will have two contributions to 
the interlayer conductivity and one is B- independent: 

a± (B) = a c ° h (B) + a^ coh {B = 0). (16) 

a coh (B) decreases with increasing magnetic field |n], |2fj| | 



work by Ho and Schofield concerning a small polaron 
model for interlayer transport. We thank them for send- 
ing us a copy of their preprint prior to submission. 



a c ° {B) 



af h (B = 0) 
y/1 + (evpcBTf 



(17) 



where vf is the Fermi velocity. If we increase B, the co- 
herent part decreases, and, therefore, T" lax would shift 
to lower values. A separation of the conductivity in two 
parts, as in Eq. (Jl6|) , has been proposed previously on a 
phenomcnological basis, in order to describe the magne- 
toresistance of Sr2Ru04 ||. (Except there a weak field 
dependence is associated with the incoherent contribu- 
tion due to Zeeman splitting). 

We have shown that a small polaron model for trans- 
port in layered systems shows a crossover from coherent 
to incoherent transport at different temperatures for in- 
tralayer and interlayer transport. It can be observed in 
ARPES, as well as in measurements of magnetoresistance 
and optical conductivity. It is sometimes suggested (or 
assumed) that the maximum in the interlayer resistiv- 
ity as a function of temperature occurs at a temperature 
jvmax determined by the strength of the interlayer hop- 
ping t±, either by fe B Tf ax - t ± or r(Tj llax ) - t x where 
r(T) is the temperature dependent scattering rate within 
the layers. However, we find that T™ ax can occur at a 
higher temperature, which is actually independent of t±, 
and instead closely related to T coh . 
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